Ultrafast and octave-spanning optical nonlinearities from strongly phase-mismatched 

cascaded interactions 
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Cascaded nonlinearities have attracted much interest, but ultrafast applications have been seri- 
ously hampered by the simultaneous requirements of being near phase matching and having ultrafast 
femtosecond response times. Here we show that in strongly phase-mismatched nonlinear frequency 
conversion crystals the pump pulse can experience a large and extremely broadband self-defocusing 
cascaded Kerr-like nonlinearity. The large cascaded nonlinearity is ensured through interaction with 
the largest quadratic tensor element in the crystal, and the strong phase-mismatch ensures an ul- 
trafast nonlinear response with an octave-spanning bandwidth. We verify this experimentally by 
showing few-cycle soliton compression with noncritical cascaded second-harmonic generation: En- 
ergetic 47 fs infrared pulses are compressed in a just 1-mm long bulk lithium niobate crystal to 17 
fs (under 4 optical cycles) with 80% efficiency, and upon further propagation an octave-spanning 
supercontinuum is observed. Such ultrafast cascading is expected to occur for a broad range of pump 
wavelengths spanning the near- and mid-IR using standard nonlinear crystals. 

PACS numbers: 42.65.Re, 42.65.Tg, 42.65.Ky, 42.65.Sf 



Ultrashort laser pulses with duration of a few optical 
cycles are fascinating because they are brief enough to 
resolve temporal dynamics on an atomic level like chem- 
ical reactions, molecular vibrations and electron motion. 
Additionally since they are very broadband and can be- 
come extremely intense, few-cycle pulses are used for 
coherently exciting and controlling matter on a micro- 
scopic level [Q. However, bandwidth limitations in the 
amplification process means that only multicycle ener- 
getic pulses are delivered from laser amplifiers, which ne- 
cessitates efficient external compression schemes Q: the 
bandwidth needed to support a few-cycle pulse is gener- 
ated by nonlinear processes, typically self-phase modula- 
tion (SPM) from a cubic Kerr nonlinearity. The nonlinear 
step leaves a strongly chirped and stretched pulse, and 
temporal compression is achieved by subsequent disper- 
sion compensation (e.g., by gratings or chirped mirrors). 

Access to ultrafast and extremely broadband self- 
defocusing nonlinearities is very attractive for femtosec- 
ond nonlinear optics. First, whole-beam and small-scale 
self-focusing problems encountered in Kerr-based com- 
pressors are avoided so the input pulse energy is practi- 
cally unlimited [Q. Secondly, with solitons both spectral 
broadening and temporal compression can occur in a sin- 
gle nonlinear material; solitons are stable nonlinear waves 
that exist as a balance between nonlinearity and disper- 
sion. For a self-defocusing nonlinearity soliton formation 
requires normal (positive) dispersion, which is predomi- 
nant in the near-IR where the majority of lasers operate. 

Self-defocusing nonlinearities can be achieved with cas- 
caded harmonic generation. This occurs when a fre- 
quency conversion process, like second-harmonic genera- 
tion (SHG), is incomplete due to absence of phase match- 
ing A/c 7^ 0, and instead the pump experiences a Kerr- 



like nonlinearity n^g^ 



(X 



-dgff/A/c, controllable in both 



magnitude and sign through Ak [Q. Cascaded SHG has 
been used for high-energy femtosecond pulse compres- 
sion |[], []-[], but has not found significant use to date. 
One reason is that ultrafast cascading is limited by the 
simultaneous requirements of being near phase matching 
and having ultrafast femtosecond response times, which 
limits the usable input wavelengths and nonlinear crys- 
tals iQ, []-Q . Therefore only multicycle compressed pulses 
were observed in /3-barium borate (BBO) crystals at 800 
nm [[], Q and at longer wavelengths (1260 nm) few-cycle 
compression was observed only in a narrow operation 
regime [Q. 

Here we show a surprising solution that has always 
been there but has been overlooked for good reason: 
when exploiting interaction with the largest x^^'' tensor 
components, strong and octave- spanning self-defocusing 
cascaded nonlinearities can be obtained even with huge 
phase mismatch and group- velocity mismatch (GVM). 
We exploit values of Ak 500 mm~^ and di2 = —500 
fs/mm that are an order of magnitude larger than pre- 
vious experiments. Ordinarily, this type of interaction is 
accessed only by quasiphase matching (QPM), providing 
an indication of the surprisingly large Ak. We demon- 
strate experimentally the utility of this process by com- 
pression of energetic 47- fs infrared pulses to 17 fs (4 cycles 
of the electric field) through excitation of self-defocusing 
temporal solitons. This occurs in an only 1 mm long or- 
dinary lithium niobate crystal. Upon further propagation 
an octave-spanning supercontinuum is observed. 

In cascaded SHG the frequency conversion from the 
fundamental wave (FW, uJi) to its second harmonic (SH, 
UJ2 = 2co'i) is not phase matched: after a coherence length 
7r/\Ak\ only weak up-conversion to the SH occurs, fol- 
lowed by the reverse process of back-conversion after an- 
other coherence length. On further propagation this cas- 
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cade process is cyclically repeated and the FW effectively 
experiences a Kerr-like nonlinear refractive index change 
An = ^casc^ proportional to its intensity /. In the strong 
phase-mismatch limit n^^^g^ —2ujid^^/c^eQn\n2/^k [Q, 
where A/c = /c2 — 2/ci, kj = njUOj/c are the wavenum- 
bers and Uj the linear refractive indices. Thus, A/c > 
gives a negative (self-defocusing) cascaded nonlinear- 
ity. However, cubic nonlinearities in the material, stem- 
ming from the electronic Kerr effect, are usually posi- 
tive > 0, i.e. self- focusing), and will therefore 
compete with the induced self-defocusing cascading non- 
linearity: to obtain a negative effective nonlinear refrac- 
tive index nl^ = n^^^^ + ^Kerr,ei we must have \nl^J > 
^Kerr el* date. Critical cascaded SHG was used for en- 
ergetic pulse compression experiments [Q []-[]. Despite 
not using the largest quadratic tensor components, in 
critical SHG the FW and SH have orthogonal polariza- 
tions and nlf^ < can be achieved by angle tuning the 
crystal close to the phase-matching point [Q. Consider 
now the normalized frequency response of the cascaded 
nonlinearity i?casc(^) = Ak/{^k^^^n'^ - di2^ + Ak) [Q, 
where di2 is the GVM and /c2^^ = is the SH group- 
velocity dispersion (GVD). Here for simplicity only up to 
second-order dispersion is included, making it possible to 
derive a threshold where incase becomes resonant [Q 

Akr = dl^/2kf^ (1) 

For > the nonresonant regime is Ak > Ak^. 
When Ak is significantly above this threshold the band- 
width Acasc = |8(A/c — A/cr)/^ 1 2 can span over an oc- 
tave. The equivalent temporal response tease = 1 /Acasc is 
subcycle making the cascading practically instantaneous. 
However, if Ak is reduced to increase the cascaded non- 
linear strength, the cascaded nonlinearity narrows spec- 
trally {Ak Akr) and then eventually becomes resonant 
{Ak < A/cr), making only weak pulse compression pos- 
sible before pulse distortion sets in [], Q. Note that the 
tunability of the cascading allows for generating a self- 
focusing nonlinearity as well, provided that Ak < 0. In 
this case the requirements for an octave-spanning band- 
width are anomalous dispersion of the harmonic (here 
< 0) and Ak < Ak^. 
Interaction through the largest quadratic tensor com- 
ponent promises an increased cascaded nonlinearity be- 
cause n^g^gc (X —dl^/Ak. This makes noncritical cas- 
caded SHG, in which the FW and SH are polarized 
along the crystal axes, very attractive. Unfortunately 
the cascading strength is at the same time reduced due 
to a complete lack of phase matching. QPM can in- 
crease n^g^g^ [|~~| by reducing the residual phase mismatch 
A/cqpm = Ak — 27r/A, where A is the QPM domain 
length, but since QPM also reduces deff (with at least 
2/7r) the breakeven value for achieving the same n^^^g^ 
as without QPM is A/cqpm(BE) = 4A/c/7r^. Let us in- 
vestigate the consequences for lithium niobate (LiNbOa, 
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FIG. 1. (Color) Performance of noncritical {0 — 7t/2) cas- 
caded SHG in LN. (a) Wavelength dependence of Ak from 
material-dispersion, the resonant threshold Akr from Eq. ([l|), 
and the breakeven QPM value AA;qpm(BE). The Sellmeier 

equation from Ref. || |] was used at T = 293 K. (b) Predicted 

\niasc\ using deff = 0^33 = —25.0 pm/V at A = 1.064 jj,m [ || , 
and using Miller's scaling to calculate the quadratic nonlin- 
earity at other wavelengths, and the predicted nKerr,ei [""l l 
using K = 3100 eV^/^cm/GW Q. (c) i?casc at A = 1.3 /xm 
when using material-dispersion, and when using QPM to get 
breakeven as well as twice as large |ncasc|- Note: incase (^) = 
Ak/[k2{^ + CJ2) — k[^''Q, — 3ki] includes full SH dispersion. 



LN): Fig. [l|(a) shows the material Ak vs. pump wave- 
length. When using QPM to reduce its large value, it al- 
ways comes with the price of having A/cqpm(BE) < Ak^. 
This is a consequence of a high Ak^ threshold caused by 
a large GVM that is typical of noncritical SHG. In Fig. 

c) the cascaded frequency response incase is plotted for 
Ai = 1.3 /im: when using QPM the breakeven case has 
a resonant cascading response as A/cqpm(BE) < A/cr, 
implying that the cascading can no longer support ultra- 
fast interaction. These properties deteriorate when QPM 
is forced to achieve a stronger n^^^^ than the material- 
dispersion case. This might have implications for ultra- 
fast cascading using QPM, such as in Refs. [| | . 

Since QPM seems generally unsuitable for increasing 
the nonlinearity, we have to use the large Ak provided by 
the material. Contrary to conventional wisdom we now 
show that in this case the cascading properties are favor- 
able and even advantageous. Firstly, as Fig. |](c) shows 
the material-dispersion case has a cascading frequency 
response that is nonresonant and extremely broadband 
(bandwidth Acasc spans 1.5 octaves). This is because 
Ak > Akr, and therefore the large GVM associated with 
noncritical interaction does not cause problems because 
Ak is so big. This holds in the entire regime supporting 
self-defocusing solitons (up to the zero-dispersion wave- 
length 1.92 /im), see Fig. [^(a). Consequently, the noncrit- 
ical cascaded nonlinearity is always nonresonant and ex- 
tremely broadband: we checked that its bandwidth spans 
over an octave in the entire plotted range and conse- 
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FIG. 2. (Color) (a) Spectra after a 1 mm LN crystal for various 
input intensities, (b) Spectra after 5 and 10 mm propagation 
using /in = 1.0 TW/cm^ (averaged over 1,100 shots). The 
spectra in (b) are shifted 20 dB with respect to each other for 
clarity. 



quently the cascading is practically instantaneous (cas- 
caded response time tease < 1 fs). Similar very favor- 
able properties were found in all other nonlinear crys- 
tals we checked as well. The question is now whether 
is strong enough when the large material Ak is 



used: Fig. 0(b) shows |n^g^g^| whose strength increases 
with wavelength due to a reduced Ak. For Ai > 0.96 fim 
it is larger than n^^^^ as calculated with the two-band 
model Q. Thus, despite the known limitations of the 
latter model, in this range an effective self-defocusing 
nonlinearity n^^^ < can be expected. 

We conducted an experiment to illustrate the ultra- 
fast and broadband nature of the cascading by using self- 
defocusing solitons to compress multicycle near-IR pulses 
towards few-cycle duration. A commercial Ti: Sapphire 
regenerative amplifier and an optical parametric ampli- 
fier generated Gaussian shaped pulses at A = 1.300 /im 
with ATfwhm = 47 fs, 59 nm FWHM bandwidth and 
200 /iJ pulse energy. Using a multiple-shot SHG 
frequency-resolved optical gating (FROG) device with a 
30 jam thick BBO crystal the measured group-delay dis- 
persion was —300 fs^, corresponding to a 42 fs FWHM 
transform-limited pulse. We found that good compres- 
sion relied on ensuring that the input pulses were not 
positively chirped (this trend was also observed numer- 
ically). A broad beam spot (0.46 mm FWHM) ensured 
negligible diffraction, and it was controlled by a reflective 
silver- mirror telescope system, and the peak intensity was 
adjusted with a neutral density filter. The crystal was a 1- 
mm-long, 5% MgO-doped X-cut congruent LN {0 = 7r/2, 
101 = 7r/2, FW and SH polarized along the vertical op- 
tical Z-axis, 10 X 10 mm^ aperture). The output pulses 
were monitored by either (a) multiple-shot intensity au- 
tocorrelator (AC) with a 100 /im BBO crystal, (b) FROG 
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FIG. 3. (Color) FROG measurements at /in = 1.0 TW/cm^: 
(a,b) FROG traces (amplitude) of the output pulse, (cl) The 
temporal and (d) spectral amplitudes and phases retrieved by 
the FROG algorithm. (c2) The calculated FROG intensity AC 
and the experimentally measured AC. 



or (c) optical spectrum analyzer (OS A). 

Figure ||(a) shows representative spectra measured by 
the OSA. For lower intensities, some spectral broadening 
occurs. When the intensity reaches /in = 0.5 TW/cm^, 
a weak shoulder appears on the blue side of the central 
lobe. At this point the AC traces evidenced pulse com- 
pression, which became strongest at /in = 1.0 TW/cm^ 
(0.13 mJ pulse energy). The FROG results in Fig. || re- 
veal more details at this intensity. The retrieved temporal 
intensity shows clear pulse compression, and a smaller de- 
layed pulse. A Gaussian fit to the central spike gave 16.5 
fs FWHM (below 4 optical cycles), and by comparing 
the fitted-curve area with the total area, we estimate an 
efficiency of 80 %. The flat phase in the time trace indi- 
cates that we are close to the soliton formation point. The 
spectral phase is also quite flat and indeed by calculating 
the transform limited pulse supported by the spectrum 
gave a 14 fs FWHM pulse, close to the observed value. 
The retrieved FROG spectral intensity in Fig. |^(d) is in 
excellent agreement with the OSA spectrum in Fig. |^, 
and similarly the FROG AC trace agrees well with the 
independently measured AC trace, see Fig. ||(c). 

Below 1.92 /im LN is normally dispersive (GVD coeffi- 
cient k^^ = +162 fs^/mm at 1.3 /im), so the combination 
of spectral broadening and temporal compression well be- 
low the transform limit of the input pulse proves that we 
have excited a self-defocusing soliton. In fact, if the non- 
linearity were self-focusing, no pulse compression could 
be achieved as such a nonlinearity would require anoma- 
lous GVD compensation after the crystal. We checked 
with an IR camera that the FW showed no diffraction. 
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FIG. 4. (Color) Numerical simulation of 2 mm propagation in 
LN with /in = 1.0 TW/cm^ The FW on- axis evolution along 
z is shown as (a) time intensity and (b) spectral intensity. 
'N' and 'A' indicate the normal and anomalous dispersion 
regimes, and the color scale is saturated for clarity. The top 
cuts show the intensities at input, z = 1 mm (the position 
indicated with a dashed white line) and z = 2 mm. (c) and 
(d) show the transverse variation along the x direction of the 
time and spectral contents at z = 1 mm and fixing y = 0. 



which can be attributed to the short crystal and broad 
spot size. Soliton formation in such a short crystal is a 
consequence of the huge FW GVD associated with using 
the largest nonlinear tensor components, which implies 
that only a short propagation distance is needed to com- 
pensate for the nonlinear SPM-induced chirp. 

Tuning the intensity around the value for optimal com- 
pression, /in = 1.0 TW/cm^, the experiment showed sim- 
ilar pulse compression and spectral traces in the range 
/in = 0.8 - 1.3 TW/cm^. Beyond 1.3 TW/cm^ substan- 
tially more complex autocorrelation traces were observed. 
For smaller intensities the temporal compression gradu- 
ally disappeared. This is typical of soliton compression 
and is consistent with numerical simulations. 

With further propagation, the spectrum continues to 
broaden rapidly; the physics behind this is similar to 
standard supercontinuum generation [Q, namely soli- 
ton fission after the initial compression point followed 
by various nonlinear interaction steps that broaden the 
spectrum. The generation of an octave-spanning super- 
continuum spectrum [Fig. ^(b)] in longer crystals is direct 
evidence of the ultrabroad bandwidth of the cascading 
process we investigate here. 

We performed numerical simulations using the 3+ ID 
slowly evolving wave approximation Q, including the 
Raman effect modeled as a 4-mode complex Lorentzian in 
frequency domain (for Raman spectrum parameters, see 
]l9[). Based on Raman peak-gain [Q and Z-scan mea- 
surements jl^, we made a detailed analysis (which will 
be published elsewhere) and found that at A = 1.3 /im the 
experimental results Q support the values n^^^^^ = 
45 X 10-20 mVW and fn = 0.50 (both ±10%), giv- 
ing n^err,el = ' fnXerr = 22.5 X lO-^^ ^2/^. In 

Fig. we used the same input pulse parameters as the 



/in = 1.0 TW/cm^ experiment. The on- axis spectrum (b) 
shows a strong central lobe, a dominant blue shoulder 
and a weaker red shoulder followed by a broad plateau. 
The on-axis time trace (a) also shows a strong spike and 
a weaker trailing pulse. The main soliton is shorter than 
what we observed experimentally, but the simulation in- 
dicates that the trailing pulse we saw is a weaker soli- 
ton formed by soliton fission (the intensity corresponds 
to an effective soliton order of around 6 for this simula- 
tion). We found that this kind of early-stage soliton fis- 
sion was mainly caused by the rather large soliton order 
combined with strong Raman effects. This fission pro- 
cess continues so after 2 mm three solitons exist, and 
upon further propagation an extremely wide supercon- 
tinuum is formed. These simulations also confirmed that 
no diffraction occurs, as was observed in the experiment. 
Moreover, concerning the spatio-temporal distribution of 
the compressed soliton, the space-time cross-section in 
Fig. 11(c) shows that it is quite homogeneous (similar to 
the Gaussian case of Ref. 0). Finally, note the optical 
Cherenkov waves that emerge in Fig. || as broadband 
waves centered around A = 3.0 /im; such waves could 
be an efficient source of mid-IR few-cycle radiation [ [ | . 
The simulations indicate the possibility of exciting self- 
defocusing solitons at shorter wavelengths as well. At 
some point, however, the effective nonlinear ity will be- 
come focusing, see Fig. 0(b). At longer wavelengths the 
Kerr, and in particular the Raman, nonlinearities will be 
less dominating, making cascading more favorable. 

Concluding, strongly phase-mismatched cascading ex- 
ploiting the largest quadratic nonlinear tensor element 
can create ultrafast and octave-spanning self-defocusing 
cascaded nonlinearities. We confirmed this by the first ex- 
perimental observation of few-cycle self-defocusing soli- 
ton pulse compression in noncritical cascaded SHG, 
which occurred in a short (1 mm) bulk lithium niobate 
crystal. Upon further propagation an octave-spanning su- 
percontinuum was generated. The noncritical interaction 
gives short interaction lengths and zero spatial walk-off, 
which represent significant improvements compared to 
previous experiments [Q []-[]. Because the nonlinearity 
is self-defocusing the pulse energy is readily scalable in 
large-aperture crystals to multi-millijoules. As it is also 
compact, cheap and works already for few- micro joule 
pulses, this compressor may replace or fill out gaps of 
standard hollow- fiber compressors |Q. 

Paradoxically, despite the strong phase mismatch a 
strong self-defocusing cascaded nonlinearity is achieved 
even without QPM owing to a huge quadratic nonlin- 
earity. However, it is exactly the strong phase mismatch 
that ensures a nonresonant octave-spanning cascaded 
nonlinearity, which is practically instantaneous, so even 
with a large group- velocity mismatch ultrafast femtosec- 
ond interaction is possible. This discovery may open a 
new era of ultrafast cascading based on strongly phase- 
mismatched interaction exploiting the largest quadratic 
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nonlinear it ies, in particular with semiconductor materi- 
als that are known for their huge quadratic nonlineari- 
ties that cannot be phase matched. It also implies that 
QPM is not needed to increase the cascaded nonlinearity. 



Thut^, rioiiit^ of Llit^ Lediiiological limlLaLioiiri of QPM ai't^ 
removed, like complexlLy in producLion, pohng quaULy, 
limited beam aperture etc. Not only that, but avoiding 
QPM is often crucial for ultrafast purposes as the cas- 
caded nonlinearity ran become resonant with QPM. 
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